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Relation between fundamental estimation limit and stability in linear quantum 
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From the noncommutative nature of quantum mechanics, estimation of canonical observables 
q and p is essentially restricted in its performance by the Heisenberg uncertainty relation, 
(A<j 2 )(Ap 2 ) > h 2 /4. This fundamental lower-bound may become bigger when taking the structure 
and quality of a specific measurement apparatus into account. In this paper, we consider a particle 
subjected to a linear dynamics that is continuously monitored with efficiency n £ (0, 1]. It is then 
clarified that the above Heisenberg uncertainty relation is replaced by (Ag 2 )(Ap 2 ) > fi 2 /An if the 
monitored system is unstable, while there exists a stable quantum system for which the Heisenberg 
limit is reached. 
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In quantum mechanics, any noncommutative observ- 
ables must possess a fundamental uncertainty due to the 
absence of their joint probability distribution. For ex- 
ample, if we estimate the position and momentum oper- 
ators of a single particle, q and p, the estimation errors 
Aq and Ap satisfy the Heisenberg uncertainty relation 
(Aq 2 ){Ap 2 ) > h 2 /A. Several type of such uncertainty 
bounds have been found in quite general formulation that 
even includes effects of measurement dQ. It is clearly 
significant to perform further detailed investigation on 
fundamental estimation limit taking the structure, prop- 
erties, and quality of a specific estimator into account. 

The quantum filter [1, 0, [H, H, 0] is a particularly im- 
portant estimator, because of its potential application to 
quantum feedback control @, ©, M, El El El EI El • 
More specifically, for a continuously monitored system, 
the quantum filter generates an optimal estimate of a 
system observable, which can be fed back to control the 
system. The estimator is recursively computed using the 
Belavkin filtering equation; this completely reflects the 
structure of the monitored system. Hence, within the 
framework of quantum filtering, the estimation limit is 
determined by dynamical properties of the system, e.g., 
the stability. 

In this paper, we particularly focus on a single one- 
dimensional particle that has a quadratic potential and 
a linear interaction with a vacuum electromagnetic field, 
the latter of which is continuously measured by a homo- 
dyne detector 0, E3, El, ES Hfl HH HI H- For this 
system, the filtering equation is reduced to the famous 
Kalman filter, and eventually the estimation error can 
be evaluated explicitly. The goal of this paper is to show 
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that, irrespective of parameters of the system, there ex- 
ists a fundamental estimation limit determined by the 
dynamical stability properties of the system. In partic- 
ular, we show that a new estimation limit on q and p 
appears if the system is unstable, while there exists a 
stable quantum system for which the Heisenberg limit is 
reached. 

We use the following notation: for a matrix A = {flij), 
the symbols A T , A', and A* represent its transpose, con- 
jugate transpose, and elementwise complex conjugate of 
A, i.e., A T = (a di ), = (a*-), and A* = (4) = (A^) T , 
respectively; these rules are applied to any rectangular 
matrix including column and row vectors. Ke(A) and 
Im(^4) denote the real and imaginary part of A, respec- 
tively, i.e., (Re(A))ij = (a i3 + a* j )/2 and (lm(A))ij = 

Kj ~ ay")/ 21 - _ 

We first review the quantum filtering theory with the 
focus on a particle interacting with a field. The inter- 
action is given by a unitary operator subjected to the 
following Hudson- Parthasarathy equation [24j : 



hdU t 



1 



dt+cdBj-cUB, 



U u U = I, (1) 



where c = c\q + c-ip. The constants Ci, C2 G C are de- 
termined according to the system-field interaction. The 
quantum Wiener process B t , which is a field operator, 
satisfies the following quantum Ito rule: 



dB t dB t = 0, 



dB\ dB t 



0. dB t dB] 



hdt, 



dB\dB\ 



0. 



In addition to the interaction, the particle is trapped in 
a quadratic harmonic potential of the form 



H = l -x J Gx = ±(0 U 



<f + guqp + gi2pg + 922P 2 ), 



Shinji'Hara@ipc.i.u-tokyo. ac.jp 



where x — (q,p) T , and G = {gij) is a 2 x 2 real sym- 
metric matrix. In the Heisenberg picture, the time- 
evolved position and momentum operators qt = U^qUt 



2 



and pt = U}pUt satisfy the following quantum stochastic 

(2) 



differential equation: 

dx t = Ax t dt + iE[CdB} - C*dB t ] 
where xt = (qt,Pt) T ■ Here, we have defined 



A := E[G + Im(G*G T )], C 



1 

-1 



Next, we consider to measure a field observable after the 
interaction. In the homodyne detection scheme, the ob- 
servable to be measured is given by 

Yl = U}(e-^B t + f**B\)U t + K(B' t + B?), 

where B' t is a noise uncorrelated from B t , and k > 
represents the strength of B' t . Also, 4> G [0, 2ir) denotes 
a phase-shift parameter that should be optimized. Re- 
defining the normalized output Y t satisfying dYf = dt, 
we have 

dY t = 2^jCjx t dt + ^(e-^dBt + e^dBj) 

+ ^T^ ) (dB' t +dB'J), (3) 

where G r := Re(e-^G) and 77 := (1 + k 2 )- 1 € (0,1]. 
Remarkably, Y t satisfies the self-nondemolition property 
[K^Yf] = 0, Vs, t for a fixed <fi, which indicates that the 
observation y t = vN{Y s | < s < t} constructs a classi- 
cal stochastic process. Furthermore, Yt satisfies the non- 
demolition condition [i^, &] = 0, [Ks,_p t ] = 0, Vs < t for 
a fixed 4>. These two properties allow us to define the 
quantum conditional expectations n t (q) — P(<7t|3^t) and 
ft tip) = P(p* \yt), which are the best estimates of qt and 
p t in the sense of the least mean square error. Follow- 
ing the quantum filtering theory, we obtain a recursive 
equation to calculate TVt(q) and fttip)'- 

c?7r t (i) = Air t (x)dt 

+ v^[|vtC r + S T Ci] [dY t - 2^jC'Jntix)dt] , (4) 

where Q := Im(e-^G) and n^x) := {w t (q), Mp)) T ■ 
Here, Vt is the symmetrized covariance matrix given by 



V t := 
Pt := 



\Pt\y t 



\(Aq t Ap t + AptAqt 



\{Aq t Apt + AptAqt) Aft 



(5) 



where Aq t := qt — i"t(<f) and Ap t :— pt — ntip) are the 
estimation errors. Vt satisfies the following Riccati dif- 
ferential equation: 



where 



V t = A'V t + V t A' T + D- ^Lv t C r CjVt : 



A' := E [G + C v Cj + (2?/ - l)GCj] , 
D := frE T [G r G r T + (1 - rfiCiCj] E. 



As Eq. ([6]) is deterministic, the quantum conditional 
expectation V t — P( Pt \ yt ) is replaced by the simple 
expectation V t — (Pt) := Tr [(p ® $)-Pt], where p is a 
system state and $ is the field vacuum state. The set of 
equations ^ and called the quantum Kalman filter 
computes the best estimate of qt and pt recursively. 

We here provide an important fact: Unlike the clas- 
sical case where the error covariance matrix is simply a 
nonnegative matrix, the canonical commutation relation 
[q,p] = ih imposes Vt to satisfy the condition 

V t + ^ E > 0, 
that yields the Heisenberg uncertainty relation 



det(Vt) > 



(Aq 2 t)(Af t )> 



(7) 



This inequality does hold regardless of a measurement 
setup. Hence the following natural question arises. 
Can the Heisenberg limit fi, 2 /4 be reached in the linear 
filtering scheme discussed above? To answer this 
important question needs a detailed investigation of 
Voo, a unique steady solution of the algebraic Riccati 
equation Voo — in Eq. ©. (If the Riccati equation 
does not have such a solution, it implies that the 
estimation fails; we do not take this bad scenario into 
account.) In particular, we aim to get a fundamental 
lower bound of det (Voo) that does not include G, G, and 
</>, because these terms completely depend on a system 
under consideration. We then obtain the following result. 

Theorem. Suppose Eq. ^ has a unique steady so- 
lution Voo. Then, the estimation error det (Voo) has the 
following achievable bounds for any G, G, and <fi: 

det (Voo) > (if G r T EG; < 0), 

dct(Ko) > ^ (if G r T EG > 0). 

Proof. The proof is done by a straightforward calcu- 
lation. Without loss of generality, we can assume that 
G r is normalized: CjC r — 1. Let G r be a unit real vec- 
tor orthogonal to G r , i.e., CjC t = 1 and CjC r = 0, and 
define 

vi := CjVooC;, v 2 := CjV^a, v 3 := CjV^a. 
Then, as (C r , G r ) is a 2 x 2 orthogonal matrix, we have 



det (Voo) = det 



CJ 
CJ 



Voo (G r ,G r ) 



ViV 3 



(6) Furthermore, let us define 



ai a 2 

a 3 a 4 

di d 2 

d 2 d 3 



CjA'C r CjA'Cj 
H T A'C C J A' c 

CjDC T CjDCj 
CjDC r CjDC, 



3 



Note that D T = D. With the above notations, the alge- 
braic Riccati equation = is reduced to 



4ry 2 

2a\V\ + 2a 2 v 2 + d\ — v 1 = 0, 



(8) 
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a 3 vi + (ai + a 4 )v 2 + a 2 v 3 + d 2 ~fr VlV2 = ^' 00 



4?7 2 

2a 3 f 2 + 2a 4 v 3 + "3 — ^-"2 = 



(10) 



Then, adding vjx^fy, -2v x v 2 x {9}, and v 2 x([10]), we 
readily obtain 

2(vi«3 - v%)(a 2 v 2 - a A v\) = d 3 vf - 2d 2 v\v 2 + d\v\. 

This together with Eq. (JSJ) leads to 



det(Voo) = -7- 



d 3 v\ — 2d 2 viv 2 + d\v\ 
Arj v\ — h(ai + a A )vi/2ri — hdi/4r) 



Note that the denominator is strictly positive from 
the assumption that the Riccati equation has a unique 
steady solution. Now, calculating d i} e.g., di = h(l — 
i])(CjY,Ci) 2 , the numerator of det(y oo ) is evaluated as 

d 3 v\ — 2d 2 v\v 2 + d\v\ 

= hv\ + h(l - r;)[(C r T EC i ) Wl - (CjZd)v 2 ] 2 > hoi 

from which we have 

ft 2 " 2 
detiVoc) > — 



v\ — h(ai + a^)v\l2r] — hd\jAr\ 

The right-hand side of the above inequality is further 
evaluated as follows. First, if a\ +a 4 = 2{r\ — l)CjT,d = 
0, which implies d\ = 0, we immediately obtain 
det(Voo) > ft 2 /4?7. Second, if CjY,d < 0, which implies 
ai + a 4 > and di > 0, we have det(Voo) > h 2 /4i]. Fi- 
nally, let us consider the case of CjEC; > that leads to 
a\ +a 4 < and d% > 0; a simple calculation clarifies that 
the function f(v) = v 2 / (v 2 + av — b), (a > 0, b > 0) satis- 
fies f(v) > 4b/ (Ab + a 2 ) when v > and v 2 + av - b > 0. 
This lower bound becomes r\ in our problem where 
a = —h(ai + a 4 )/2?y and b — hdi/Ar/. As a result, we 
obtain det(Voo) > fi, 2 /4 in this case. The achievability 
of the above lower bounds is discussed in the example 
part. ■ 

We now give a physical interpretation to the sign of 
CjY,Ci. To do this, let us focus on the matrix A, which 
corresponds to the drift term of the quantum dynamics 
@ and the filter ([J}. The characteristic polynomial of A 
is A 2 + 2(C I T EC i )A + (C r T ECi) 2 + det(G) = 0. Hence, A 
has two stable eigenvalues if and only if the conditions 

C r T £Q > 0, (C r T SCi) 2 + dct(G) > (11) 

are satisfied. The latter condition is easily attained by 
making the coefficient of C (i.e., the interaction strength) 
sufficiently large, if the former condition is already sat- 
isfied. Therefore, under the condition C^SCi > 0, both 



the quantum dynamics and the filter are (asymptotically) 
stable in the sense that, roughly speaking, those trajec- 
tories are constrained around x = and 7r t (i) = 0. This 
implies that the fundamental estimation limit h 2 /Aij can 
be violated if the dynamics we aim to track is stable. 
Combining the theorem with the above discussion, we 
deduce the following fact: 

ft 2 

(Aql Q )(Apl c ) > — (if the system is unstable), 



> 
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(if the system is stable). 



Remark 1. In practice we cannot construct a perfect 
measurement apparatus with 77 = 1. Thus, the condition 
CjUCi > is clearly preferable from the estimation 
performance viewpoint. Actually, for example when 
C^SCi < and 77 = 1/4, the estimation error is lower 
bounded by ti 2 , i.e., (Aq 2 < ,)(Ap 2 x ,) > h 2 , which is 
much bigger than the Heisenberg limit ft 2 /4. However, 
the sign of C^SC; cannot be changed by tuning the 
Hamiltonian matrix G and the phase-shift <j>. (Note that 
C r T £Ci = Re(C) T SIm(C).) In other words, only the 
interaction term C is the crucial factor that determines 
the estimation limit. 

Remark 2. The Hamiltonian of the form H = 
x T Gx/2 — x T 'SBut, where B £ M. 2 , allows that the system 
dynamics 

dx t = Ax t dt + Bu t + iS[CdBl - C*dB t ] 

can be controlled using a feedback input u t £ 3^t- For 
example, the quantum linear quadratic gaussian (LQG) 
controller effectively stabilizes the system. However, 
any control input cannot reduce the estimation limit, 
because the error covariance matrix Vt obeys the same 
Riccati equation ^ without respect to B and Ut- 

We will show that the two bounds in the theorem are 
tight in a sense that there exists at least one example 
where the equality holds in each case. 

Example 1. Doherty et. al. considered in [l|| a sin- 
gle particle system with the following harmonic oscillator 
potential and the interaction with strength a > 0: 



H = q 

2 q 

This corresponds to 
■ 2 



2m 



G = 

C r = 



raw 




1/ra 

2a cos 4> 




C = - 
Ci = 



'2a 



2aq. 



1 



'2a sin ( 




First, we remark that the Ricatti equation © has a 
unique steady solution Voo for all the parameters. Then, 



4 



due to C^SCi = 0, the estimation error is bounded by 



h 2 

det(Ko) > — 
477 



(A^)(A^> > -. 



Actually, the drift matrix A has eigenvalues ±iu>, imply- 
ing that the particle is oscillating with frequency u), and 
thus that the system is not stable. Furthermore, in this 
case, we can obtain a simple explicit form of det(V 00 ): 



det^) 



\COS z / 



4r/ Vcos 



which attains h 2 /4f] when <fi — 0. Therefore, the lower 
bound ti 2 /4r] is indeed achievable. In particular, when 
6 = we have 



fi 2 /4?7 



4r/ yj r \ + n + h 2 / T 4r] 

where r% = hmu> 2 /8i]a. Thus, in the limit of r± — > oo 
the estimation error satisfies the minimum uncertainty 
relation (Aq 2 <} )(Ap 2 a ) — h 2 /Arj : which further attains 
the Heisenberg limit {Aq 2 ){Ap 2 xl ) = h 2 /A only when 
r) = 1. 

Example 2. Wiseman and Doherty considered in [l9[ 



an atomic system in a damped cavity containing an on- 
threshold parametric down converter that realizes 

j3 

H=-{qp+pq), c = j(q + ip), 
where (3 > and 7 > are parameters. We then have 



Cr = 7 



a = 7 



The Ricatti equation ([6]) has a unique steady solution 
under the condition f3+j 2 > 0, which is already satisfied. 
Then, due to CjSC; = j 2 > 0, the estimation error is 
lower bounded by ft 2 /4 from the theorem. This bound 
is achievable as in the former example. Actually, when 
4> = 0, the off-diagonal term of Voo is zero, and eventually 
we have 



(aO(a/ c 



tf_ \Jr\ + 2(277 - 1>2 + 1 + r 2 + 2?7 - 1 

8r?' l+r 2 



where r2 := (3/ J 2 - Hence, when 7 — ► 00, which implies 
that the interaction strength is very large, the Heisenberg 
limit h 2 /4 is reached. 
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